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Abstract 

The effects of three-dimensional perturbations in two-dimensional turbulence 
are investigated, through a conformal field theory approach. We compute scal- 
ing exponents for the energy spectra of enstrophy and energy cascades, in a 
strong coupling limit, and compare them to the values found in recent exper- 
iments. The extension of unperturbed conformal turbulence to the present 
situation is performed by means of a simple physical picture in which the 
existence of small scale random forces is closely related to deviations of the 

exact two-dimensional fluid motion. 
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I. INTRODUCTION 



It has been recognized that turbulence, with its manifold experimental realizations, is 
one of the challenging problems to which very different approaches, ranging from pure math- 
ematics to engineering applications, have been developed in an interesting complementary 
way. One of the most important methods to study turbulence is, in fact, the formulation 
via field theory, based on its relationship with stochastic partial differential equations ]!],[§. 
However, such technique is far from being well established and complete, so that new ideas 
and important improvements are constantly appearing in the subject. 

Recently, Polyakov suggested that non-unitary minimal models of conformal field theory 
could be used to describe the physics of two-dimensional turbulence |§ . The advantage of 
this proposal is that one can deal in a controllable way with a set of anomalous dimensions 
and short distance products. An infinite number of inertial range exponents follows from 
this approach and one of the still open problems is how to find "selection rules," which 
would define the experimentally relevant minimal models or the connection between them 
and statistical characterizations of the random forces acting on the system. These ideas 
have attracted the attention of many authors, and generalizations have been investigated, 
as, for instance, possible boundary effects 0, alternative physical pictures for the enstrophy 
and energy cascades |J, and magnetohydrodynamic turbulence ||. 

We will consider, in this paper, the problem of conformal turbulence including in its 
formalism the influence of three-dimensional effects. Our motivation comes from a number of 
experimental studies, in which approximately two-dimensional fluids were observed, together 
with the unavoidable presence of three-dimensional perturbations [p~0| P~3f| . It was verified 
that a quasi two-dimensional fluid is perturbed by small scale forces originated from the 
the degrees of freedom related to the direction perpendicular to the plane of motion. We 
will take this fact into account, noticing that there are also compressibility effects which 
cannot be neglected in an effective two-dimensional theory of the perturbed system. A 
generalization of the conformal approach will be devised and new inertial range exponents 
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will be obtained here, in reasonable agreement with the experimental data. 

This paper is organized as follows. In the next section we briefly review the most im- 
portant and practical aspects of conformal turbulence, in order to make the paper as self 
contained as possible. In section III, we discuss some of the experiments carried out to inves- 
tigate two-dimensional turbulence. This will motivate us to define an effective (and infinite) 
set of stochastic partial differential equations which represents a quasi two-dimensional fluid 
under the influence of three-dimensional perturbations. The conformal approach is then 
introduced in order to solve the Hopf equations for the turbulence problem. Furthermore, 
the constant enstrophy and energy flux conditions are also studied. Explicit solutions are 
found and described in section IV and, in section V, the problem of boundary effects is 
discussed. Finally, in section VI we comment on our results and on possible directions for 
future investigations. 



II. CONFORMAL TURBULENCE 



The minimal models of conformal field theory |14]] are generically defined by a pair of 
relatively prime numbers, (p, q), with p < q. These models contain a subset of (p— 1)/2 
scalar primary operators, ^(m.n); labelled by 1 < m < p and 1 < n < (q — l)/2, if p 
is even, or 1 < m < (p — l)/2 and 1 < n < q, otherwise, having dimensions A( mn ) = 
({pn — qm) 2 — (p — q) 2 ) /4pq. The reason for the choice of scalar operators is that we will be 
dealing with isotropic correlation functions in the turbulence problem. The operator product 
expansion (OPE) of two primary operators ^(n.si)^) and i>(r2,s 2 )( z ')i with \z — z'| — ^ is 
written as 

VWotolW*') = E (aa)( A <-3)-A(n,n)-A ( , 2 , S2) ) £ Cfc^y^^y 

(73,53) (n,m) 

L_ ni ...L_ nfc L_ mi ...L_ m; a^ n a^ m T/>( r . 3iS3 )(2:) , (2.1) 



x 



where \ri — r 2 | + 1 < r 3 < min(ri + r 2 — 1, 2p — r\ — r 2 — 1), \s\ — s%\ + 1 < S3 < 
min(si + S2 — 1, 2q — s\ — s 2 — 1) and we have introduced, in (|2.1|) , the Virasoro generators of 



conformal transformations, L_ n and The interest in these models is related not only 

to their finite number of primary operators, but also to the fact that their dimensions and 
the form of short distance products are completely known. 

Let us look now at the problem of turbulence in two dimensions and show how it may 
be matched [[J with the above operator structures. The motion of an incompressible fluid is 
assumed, even in the turbulent regime, to be described by the Navier-Stokes equations for 
the velocity field, 

d t v a + |tf a7 - vpdpv*, = f a + vd 2 v a , (2.2) 

where f a represents a random force acting at large scales, determined by a characteristic 
length L, and v — > is the viscosity, associated to the small scale where dissipation effects 
come into play, yielding a natural UV cutoff for the system. In terms of the stream function, 
ip, related to the velocity field by v a = epadpip, we may write the following equation for the 
vorticity field, u = d 2 ip, 

d t uj + e a/3 d a ipd 2 dpip = e a pd a fp + ud 2 u . (2.3) 

One of the fundamental problems of turbulence theory is to find solutions of the Hopf 
equations, for statistical averages over realizations of the velocity field, 

d t [< u(x u t)u(x 2 ,t)...u(x n} t) >} = , (2.4) 

where the time derivative is expressed through the use of equations (|2.3|) . In the inertial 
range, the standard view of the problem is that both forcing and viscosity terms may be 
neglected in order to formulate an effective set of Hopf equations. Considering, further- 
more, the convection term in ( |2.3j ) as a vanishing point-splitted product of fields, that is, 
^ z _ z/ \ = ^(dz' /a)e a /3d a ip(z)d 2 df3ip(z') — > 0, when \z — z'\ — ► 0, we would have, then, an exact 
solution of ( |2.4| ). A concrete realization of this possibility may be achieved if we regard the 
stream function if} as a primary operator of some conformal minimal model. In this case we 
may use all the available information on operator dimensions and OPE's to extract physical 
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results from the analysis of the problem. According to this assumption, let be the primary 
operator which has the lowest dimension, A<p, appearing in the OPE ipip, between fields 
with the same dimension A-0. Taking a = \a\ exp(id), we will have, thus, 

lim / —e a/3 d a 'ip(z)d 2 dpip(z') 

\a\->0J\z-z'\ = \a\ a 



m Mz,z) 



J d9 [dld a d z - d 2 a d- a d z ] (aa)( A ^ 2A ^^C {n;m} L_ ni ...L_ nfc L_ mi ...X_ m! aS"aS' 

L_ 2 L\ - L_ 2 L 2 A , (2.5) 



as the dominant contribution in this short distance product. It is important to note that 
in order to get ( |2.5| ) it was necessary to set C{i ;2 } = C{2 ; i} and Cj^^m = Cm^iy, as it 



follows from the pseudoscalar nature of the e factor above. We see, then, that ( |2.5| ) vanishes 
with \a\ — > if 

A0 > 2A^ , (2.6) 

which is one of the constraints that the chosen minimal model has to satisfy. An additional 
constraint comes from the condition of a constant enstrophy or energy flux through the 
inertial range, meaning that < ii)(x)u)(0) >~ a; or < v a {x)v a (0) >~ x°, respectively. In the 
case of a constant enstrophy flux, we have 

(A^2Av>) < ]^L_ 2 ll x - L^Ll,) 0(x)] d 2 ij;(0) > . (2.7) 



The correlation function at the RHS of (|2.7|) is now evaluated by means of a purely dimen- 
sional argument, as £ J ~ 2 ( A 'I , + A ^+ 3 ) ) which makes sense if one thinks that there is an effective 
IR cutoff in the theory at the scales where the forcing terms act. Imposing ( |2.7|) to be 
independent of L, we get 

A0 + A^ + 3 = O. (2.8) 

In the case of an energy cascade, the argument is the same and the constraint turns out 
to be A0 + Aip + 2 = 0. It is known that there is an infinite number of minimal models 
compatible with (|2.6|) and ( |2.8|) 0. The general belief, and still an open problem, is that 



there may be universality classes, associated to the statistical properties of the forcing terms, 
which would single out one or another of the possible solutions. 

An alternative analysis of conformal turbulence regards the existence of boundary effects 



on the vacuum expecation values (VEV's) of single operators in non-unitary theories ||T5 
In this case, one has to consider the OPE between 4>(x) and i/j(0) in (|2.7|) , picking up the 
most relevant operator, let us say, x- Now, ( |2.8|) is modified to A<p + Aip — Ax + 3 = 0, with 
an analogous change for the constant energy flux condition. Some of these further solutions 
(in the enstrophy cascade picture) were obtained in ref. Q|. 

The connection of the conformal approach with real experiments or numerical simulations 
is made through the computation of inertial range exponents, which describe the decrease 
of energy in the region of higher Fourier modes. In the situation where VEV's of single 
operators vanish, the inertial range exponents are given by 4Aip+l and, in the opposite case, 
by AAip— 2A0+1. A good agreement has been reached between the former possibility, for the 



the direct enstrophy cascade case, and numerical simulations pl6| , |TT|| of the two-dimensional 
Navier-Stokes equations. 

III. THREE-DIMENSIONAL EFFECTS 



In a series of interesting experiments, Hopfinger et al. ||T0|--|12| studied the turbulence 
phenomenon as it happens in a rotating tank, where at its bottom there was an oscilating 
grid responsible for perturbations of the fluid motion. According to the Taylor- Proudman 
theorem |18|-^(J a rotating fluid tends to behave as if it were two-dimensional and in fact this 
was observed in the form of coherent structures (vortices) organized in the direction parallel 
to the rotation axis of the tank. However, "defects" in the vortices were seen to propagate 
from the very turbulent region at the bottom of the tank up to the effectively two-dimensional 
system. The essential picture extracted from these observations is that the fluid should be 
best described in terms of two-dimensional equations containing not only large scale forcing 
terms but also small scale random perturbations, originated from either vortex-breakdown 
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or soliton pulses propagating along vorticity filaments. The experimental data suggested 
then the existence of an inertial range, likely to be related to a direct enstrophy cascade and 
well approximated by E(k) ~ k~ 2 ' 5 , which represents a less steep energy spectrum than the 
one obtained by Kraichnan pi| , E(k) ~ k~ 3 , or even other proposals p2| , p3| , not excluding 



conformal turbulence [Q. This puzzingly result is presently understood to be due only to the 
measurement techniques used in the experiments, based on the analysis of the dispersion of 
suspended particles in the fluid [Q. More recently, similar experiments were conduced by 



Narimousa et al. fl3| and direct measurements of the turbulent velocity field were recorded. 



The results pointed out the existence of a possible energy spectra E(k) ~ k 5//3 at lower 
wavenumbers, in agreement with the conjecture of an inverse energy cascade ||21|| , and a range 



at higher wavenumbers, where E(k) ~ k~ 5 - 5±0 - 5 . In this region, the spectral slope was seen to 
depend on the controlling external conditions, with results varying from E(k) ~ A; -5 up to 
E(k) ~ A; -6 . It is worth to note that a spectral law E(k) ~ k~ 5 follows from Rhines theory 
of /?-plane turbulence and, alternatively, is closely approximated by some solutions of 
the constant enstrophy flux condition in the conformal approach, like the minimal models 
(9,71) or (11,87). 

The variation of exponents obtained in the experiments may have a theoretical counter- 
part in the existence of a set of operator anomalous dimensions, making it interesting to 
analyze the problem from the conformal field theory point of view. It is clear, however, that 
the inertial range exponents, found in ref . [|J , cannot reproduce the experimental situation. 
We believe that the important ingredient, missing in the previous conformal approach, is 
precisely the existence of three-dimensional perturbations, which must be taken into account 
in any realistic model of a quasi two-dimensional fluid. 

In view of the above considerations, let us write the two-dimensional Navier-Stokes equa- 
tions as 

d t v a + v p dpv a = vd 2 v a + /« + gf&> - 8 a P , (3.1) 
where and ffp are stirring forces defined at large (L) and small (fi « L) scales, 



respectively. The dimensionless constant g represents, roughly, a coupling with the three- 
dimensional modes of the fluid. We assume that the dissipation scale, a, is related , in 
principle, to the other scales of the problem asa<</i<<L. This means that even though 
the perturbations act at very small scales, when compared to the macroscopic size of the 
system, they are still much larger than the scale where dissipation occurs. 

An important point here is that the condition of incompressibility, when formulated in 
three dimensions, reads diVi + c?2f2 + d^v^ = 0, suggesting that the "projection" of this 
constraint to the two-dimensional world has to be given by d a v a = O(g), in the framework 
of equations (|3.1|) . The velocity field may be described, then, by means of a stream function, 
ip, and a velocity potential, 0, as 

v a = epadpi/; + gd a (j) . (3.2) 

It is of further interest to study, besides the vorticity u, the divergence of v a , given by 
p = gd 2 (p. An exact, although infinite, chain of equations may be generated if we expand 
ip and <p in powers of g, substituing them into (|3.1| ) and collecting the coefficients of the 
obtained series. Defining, in this way, 

oo oo 

n,/,(n) ,., _ V* n n uj( n ) 



n=0 n=0 



oo oo 

E 9 n <P {n) , P = E 9 n+1 P in) , (3-3) 

n=0 n=0 



we get the following set of coupled equations, 



n— 1 



i) <W n) + E e af3 d a tp (p) dpd 2 ^ n - p) + E [d^dpd 2 ^-^ + d 2 <p {p) d 2 ip {n - p - 1] 



p=0 p=0 
,(2 
IP 



+ e aP d a ff ] 5 ritl 



n) d t u;^+e aP d a ^d fi d 2 ^ = v&u>® + e aP d a fU , 

n-1 

in) d tP {n) + E [d a d^ p) d a d^ n ' p - l) + d^dadP^-v-V 

p=0 



+ E ^Bapdpd^dM^ + e^d^^d/sd 2 
p=0 
n+1 

+ E [d a d^ {p) d a d^ {n - p+l) - d 2 ^ p) d 2 i) {n - p+1) ] = ud 2 p {n) 



p=0 



iv) d tP ^+2d a d^d a d^ + 2t a pdpd a <P^d a d a ^ + e^d^dpd 2 ^ 

-2d 2 ^d 2 ^ = vd 2 p^ + d a f® , (3.4) 

and, finally, the constraint of incompressibility for the g-independent part of the velocity 
field, which defines the pressure term, 

(cW (0) ) (<W (0) ) - 0V°W O) = d a fV - d 2 P . (3.5) 

In the above expressions, n > 1. We have obtained, therefore, a set of stochastic partial 
differential equations. In a statistical description, reflecting a stable asymptotic limit for the 
correlation functions of u and p, Hopf equations may be straightforwardly written as 

N M 

d t <Y[cu^\x u t) n P (n) (x 3 ,t)>=0 . (3.6) 

i=l j=N+l 

We observe now that in ( |3.4| ), equation ii) is identical to the one which corresponds to 
an unperturbed (g = 0) two-dimensional fluid. This means that the field ^ (0) will be 
related to an enstrophy or energy cascade, even in the presence of three-dimensional effects. 
This field plays the role of an external random variable in the other equations, since its 
dynamics is independent of the other components ip^ or to the field <p (i n general, the 
subset {ip(°\ ^ j ••■) ?p( n \ 4>(°\ 4>^\ (p( n ~^} contains fields which act like external random 
perturbations in the equations for ip^ and with p > n + 1). Considering that ( |3.4|) 

gives relatively complex equations, the analysis of the problem might seem hopeless, being 
perhaps adressed only to a numerical treatment. However, we can extend the conformal 
approach, applied previously to the unperturbed case, to find here solutions of the Hopf 
equations. Our basic assumption is that not only ip^ but also the other components in the 
power expansions of ip and are primary operators which belong to some minimal model 
in a conformal field theory. It is necessary, then, to define a scale £, possibly associated to 
intermittency effects, which allows us to write the following dimensionally correct expansion, 

oo 

n=0 

oo 

n=0 
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where ip^ and (p^ have dimensions Aip^ and Acj)( n \ respectively. 

The introduction of a scale I in (|3.7| ) means that the perturbed system exhibits a breaking 
of scale invariance in the inertial range. We will see later that this phenomenon is signaled 
by the existence of constant enstrophy or energy fluxes which depend on the small scales 
of three-dimensional perturbations. It is conceptually important to understand the physical 
origin of I. A clue for this comes from the structure of couplings between ip^ and the other 
fields, as expressed in equation ( |3.4| ). As we have already observed, ip^ is effectively an 
external field in the equations for ip( n > (with n > 1) and (n ) (for any n). In this way, it 
is plausible to have a relation between £ and the scales involved in the dynamics of tp^ 1 '. 
Now, if we consider the turbulent limit of the equations for ip(°\ corresponding to v — > 
(or, alternatively, a — > 0), we are left essentially with the correlation length L of large 
scale random forces. A simple choice, thus, is to consider i = L. In this respect, one may 
observe that the small scale /i could also be used in the definition of I. We have, however, 
physical reasons to believe that this does not happen: \x is related to the forcing terms in the 
equations for ip^ 1 ' and (p^\ which we expect to be irrelevant when compared to the nonlinear 
convection terms in the range of wavenumbers given by \k\ « 

It is interesting to note that there is an analogy between our problem and the statistical 
mechanics of second order phase transitions for a system close to its critical point. In this 
case, one can study deviations of the critical temperature T c by means of an expansion in 
(T — T c ) and through the use of the operator structure of the critical theory . Here, in 
the turbulence context, the "critical theory" is just what we get when g — > 0. 

In order to simplify the notation we will keep using (|3.3|) , with the above observations in 
mind. We are interested to obtain possible combinations of primary operators, in eq. Q3.7D , 
that would not affect, in the limit \i — *■ 0, the constant enstrophy or energy fluxes, obtained 
from the dynamics of the field Within this point of view, it is important, therefore, 

to consider short distance products of a certain number operators, as it follows from fl3.4p . 
Taking two generic primary operators, and 0% (f° r example, <fi^ and ip^'^)-, with 
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dimensions AO^ and AO^\ respectively, we may write 



(^AA < g ,p ' c ) - \< >''•'- \ <•)"•'' 

... ) ■- J \ y 

(p) 



1 W(z > f)Or ) (^^ = (aa) ,Al -' , '''" A ' ,,, ''" Aa 



x 



X^{(ni,Lin fc );( m i,...,m f )}^-«i-"^ Z ) > ( 3 ' 8 ) 



where Aq'^q is the primary operator with the lowest dimension in the above OPE. The 



short distance products appearing in ( |3.4| ) are listed below together with the conformal field 
theory representation, obtained after straightforward computations: 

.d a d^d a d^ ~ d a ^d a ^ ~ (aa)(^-^ W -^ , - 1 )L_ 1 L_ 1 ^ 
v r 2 r r 2 1 4(Pip') 

.d a dffll>V>d a dffll,W ~ d^d^ ~ (aa) K^-a^-a^-i) j (3 9) 

with > 0, except in the last relation, for the product of the type ?/>?/>, where p + p' > 1. 
From ( |3.9| ) we see clearly that Hopf equations are satisfied if 

i; A<f - 2A^ (0) > , 

z*; aaJ/ } - A^ (p) - A^ (p,) - 1 > , 
aaJ/ - av> w - A0( p ') - 1 > o , 

iv) AA { ff' ] - A0( p ) - A0( p ') - 1 > , (3.10) 

with p + p' > 1 in ii) and > in in) and wj. These equations are the first step in 
the generalization of unperturbed conformal turbulence, in order to deal with a larger set 
of primary operators. We will find more inequalities, restricting, then, up to some extent 
the number of different operators allowed in the theory. Let us now write the conditions 
for constant enstrophy or energy fluxes through the inertial range. Here we will assume 
that vacuum expectation values of primary operators are zero. The alternative possibility is 
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discussed in the next section. The case of a constant enstrophy flux requires, as commented 
before, that we compute < lj(x)lj(0) >. From relations ( |3.4| ) we get 

< cb(x)u(0) >= 

oo n 

E E 9 n+m < {e^d^dpd 2 ^^ + g 

n,mp=0 

xd 2 ^ m) (0) > . (3.11) 

In the above expression we may define a "large scale" part as the one which depends solely 
on the field ip^ and a "small scale" part, involving the fields (f)^ and the other components 
of ip. Regarding the large scale part, we already know, from the study of the unperturbed 
case, that the constant enstrophy flux condition is 

AA§f + A^ {0) + 3 = . (3.12) 

It is natural to assume, like in the case of unperturbed conformal turbulence, that the 
correlation functions in the small scale part may also be evaluated by means of a dimensional 
argument, where, instead of using the typical large scale parameter L, the correct choice 
turns out to be the small length scale ft. Assuming, furthermore, that fi — *■ leads to a 
well defined limit, we just require the powers of fx in the most relevant terms belonging 
to the small scale part of ( |3.11| ) (contributions which have the lowest power of aa) to be 
non-negative numbers. This discussion may be restated by saying that we will need to select 
one or more of the following conditions, 

i) AA% p,) + A^") + 2 < , 

ii) AA^f + A</> (p,,) + 3 < , (3.13) 

according to the analysis of the dominant terms in the small scale part of < uj(x)uj(0) >. In 
the derivation of ( |3.13| ) we have used the OPE's computed from the Hopf equations, given 
by (|3.9| ). An additional care must be taken if it happens that A^f ' = tjj( p "' ) or A^ = ip^ p "^ 
for some values of p,p' and p" . In this circumstance it is necessary to have Aip^'^ = —3/2 or 



dp^dfj d 2 i/j (n - p) + <9 2 (p) dV (n ~ p) } 
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Aip( p "^ = — 1, respectively, to assure spatially independent correlation functions and hence 
a constant enstrophy flux. 

Let us turn now to the case of a constant energy flux. We have 



oo n 

-EE 9 n+m < {g%<P {p) dpd a <P^ + g 

n,mp=0 

+d a d~ 2 [«9V 0) <9y o) - d a dp^d a d^ 
x (^W(o) +5W H(o)) > . 



5n,o + e a pe ia d^dpd^ n -^} 



(3.14) 



Here we cannot refer immediately to the Hopf equations and formulate a set of conditions, 
as we did in the constant enstrophy flux case. There are, in ( |3.14| ), OPE's which do not 
appear in (|3.4|), viz. 



«/(p)««/(pO / ,(AA<7' ) -A^)-A^')-l) 

*e a(3 e lz d (T il)^>d f5 d 1 ^ p > ~ (aa)\ w > ' L_ X A$$ 
re^e^d^dpd^ + d z d~ 2 \d 2 ^ 0) d 2 ^ 0) - d a d^d a d^ 



(a^ 0) -2AV>( > 



L_iL_ 2 - Ad-^L'^L 



A 



(0,0) 

Iplp 



.fl^W" ~ ( a)K*' > -^"- a *"' > - I ) L _ 1 ^) . (3.15) 



In ( |3.15| ) we have e 7Z = (e 7 i — ie 7 2)/2. Using the above point-splitted products, we can get 
the constant energy flux conditions. First, the equation following from the large scale part 
of < (0) >, 



A^ 0) + AV (0) +2 = 



(3.16) 



and then the inequalities which come from the small scale terms, 



i) AA$f + + 1 + W,o < o , 
n) AA^f + A0(p") + 1 + Vp'.o < , 
in) AA^f } + A^ (p,,) + 1 < , 
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iv) AA { ™' ] + A0(p") + 1 < 



v) AA^f + A0^") + 1 < 



1 (p.. 



vi) AA% P>) + A^ (p,,) + 1 < , (3.17) 



where, analogously to the enstrophy cascade case, only a subset of ( |3.17|) has to be consid 



ered. Supplementary relations, like the ones obtained after equation ( |3.13| ), are in order, to 
avoid possible ^-dependent correlation functions. We are led, here, to 

i) AV» = -1, if <f = V»W , 
t<;A^W = -l,ifi4^ 0) = ^«, 

m) Atfj^ = -1/2, if A^f ) = ^ p "\ for p + p' > 0, or Affl = 
or A™ = , 

iv) A^ = -1/2, if A§f> = for p + p' > 0, or A% p>) = ^ 

orA$f = ^'\ (3.18) 

Once we have some solution at hand, derived from the conditions obtained here, we 
may associate inertial range exponents to each one of the fields ift^' and 0^, expressed 
by AAip^ + 1 and AAcp^ + 1, respectively. From these values, we have to select the one 
which will appear effectively in experimental situations. This problem is investigated in the 
following section. 



IV. ANALYSIS OF THE CONSTANT FLUX CONDITIONS 

We have obtained so far all the conditions necessary to find minimal models related 
to an enstrophy or energy cascade in a quasi two-dimensional fluid. In order to explore 
them, the first observation we can make is that these models must belong to the infinite set 
of solutions found in the former study of unperturbed conformal turbulence. This follows 
directly from the conditions which depend only on ip(°\ A strategy of computation could 
be, thus, just a numerical analysis of all possible combinations of fields for these previously 
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known minimal models. As straight it may sound, this approach is hardly useful when the 
number of primary operators becomes large, a fact that happens already for the first few 
minimal models. 

A more interesting computational scheme is provided if we look for solutions of the form 

</> = M<h , (4-1) 

where / a (0) = /&(0) = 0, that is, we are considering solutions with ip^ = ipi, for p > 1, 
and (j)^ = 00; f° r any V- This approach is valuable since a little reflection shows that if 
it is impossible to satisfy the constant flux conditions through any pair of fields ip\ and 
00) then there are no further solutions for the model under consideration. All our task is, 
therefore, to consider the set of minimal models representing conformal turbulence without 
perturbations, from which the fields -0 O ma y be immediately obtained, and add, according 
to the new constraints associated to three-dimensional effects, the fields ip\ and 0o- 

In the study of the inertial range exponents, we may think of, at least, three limits for 
f a ,b(g): a) g -> 0, that is, f a ,b(g) -»• 0, b) f a ,b(g) ^ 1, and c) g » 1, which may be defined 
as a "strong coupling" regime. In the first case, the perturbations play a negligible role 
and everything is described by unperturbed conformal turbulence. A competition between 
exponents shows up in the second case, where the less steep spectral slope will be the 
most relevant in the limit of higher wavenumbers. We see, in this way, that cases a) and 
b) cannot give any of the steeper spectral slopes observed in real experiments. The third 
case is, in fact, where we have some hope to find a relation with experimental results. 
It would be unphysical to have f a ,b(g) 0, for large values of g, since in this limit we 
would recover the unperturbed system. Also, it is unlikely to have f a ,b(g) —> const.: taking, 
for instance, gaussian random forces /£> and f£\ with < f£\x,t)ff\x\t') >= and 

\)S(t — t'), it may be proved, from the retarded nature 
of the diffusion propagator, that < f^(x,t)v/3(x',t) >= gD a/3 (\x — x'\), yielding 
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d_ 

dg 

Let us assume, thus, that f a ,b(g) diverges as g — > oo. This means that the inertial range 
exponent derived from ip may be discarded and we have to analyze only the competition 
between the exponents obtained from ipi and O - 

We performed an investigation of the first six minimal models for both the enstrophy and 
energy cascade cases. In the enstrophy case we found solutions for all the models studied. 
They are represented in table 1 and in figure 1. In table 1, we show the fields ipi and <po for 
the minimal models (2,21), (3,25) and (3,26), together with their associated inertial range 
exponents. As the number of solutions became larger, we had to represent the other three 
models, (6,55), (7,62) and (8,67) in figure 1, where we plotted the most relevant exponents, 
found from the competition between ip 1 and O > i n the strong coupling regime. We observe, 
from the results, that there is a good agreement with experimental verifications, with the 
only considerable deviation occuring for the very small set of two solutions for the model 
(2,21). The solutions, excluding the model (2,21), were organized in table 3, where values of 
mean exponents and standard deviations are described. It is clearly seen that the perturbed 
exponents are in general lesser than the exponents of the unperturbed fluid. 

In the energy case, an interesting fact happened: most of the models studied did not yield 
any solution for the fields ipi and 0o- Only the model (10,59), represented in table 2, gave 
solutions, all of them with inertial range exponents close to —3.0, which do not support the 
conjecture of a Kolmogorov exponent, —5/3, for the range of lower wavenumbers. However, 
more theoretical and experimental work is necessary in order to arrive at a conclusive answer 
on this point. 

V. BOUNDARY PERTURBATIONS 

It is worth to understand what happens when boundary effects are supposed to have some 
influence in the problem of conformal turbulence. Below, we obtain the set of conditions 



<f^(x,t)v p (x',t)> 



D, 



X 



X'\) 



1 . 



(4.2) 
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needed to account for it, leaving a numerical analysis for future investigations. 

The basic modification here is that we have to study further OPE's in the conditions of 
constant enstrophy or energy fluxes, found in section III, since now VEV's of single operators 
do not necessarily vanish. In this way, let us define the primary operator ^[0102)03 as ^ ne 
one with the lowest dimension appearing in the OPE (0 ( f ) Oi P ' } ) 0f\ where the product 
of and was computed first. The conditions we are looking for must be obtained 
from the analysis of the x-dependence of the dominant terms in ( |3.11| ) and (|3.14j) . In the 



situation of a constant enstrophy flux, the large scale part of Q3.11| ) gives 

A^+A^)-AA;X + 3 = 0, (5.1) 

which is nothing else than the condition established in section II, in a diferent notation. On 
the other side, the small scale part of Q3.11| ) gives one or more of the following conditions: 

AA^f + A^") - AA<*/2" + 2 = 0, 



ii) AA$f + A^") - AAffiff + 3 = 0. (5.2) 



A similar analysis for the case of an energy cascade yields, for the large and small scale parts 
of fl3.14| ), respectively, 

A<f + A^°) - A<X + 2 = (5.3) 



and 



i) AA$f + A^") - AA^f + 1 + W,o = , 

ii) AA$f + A^") - AA^f + 1 + W,o = , 
Hi) AA$f + A^") - AA^f + 1 = 0, 

iv) AA^f + A<j>^ - AA^f + 1=0, 

v) AAff + A^") - AA^f + 1 = 0, 

vi) AA$f + AV>"> - AA^f + 1 = 0. (5.4) 
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The computation of inertial range exponents is also modified. We now have to consider all 
possible combinations like ip^ip^ and (jy^cp^ in the evaluation of the velocity-velocity 
correlation function. The observed inertial range exponent must be obtained from 2A^^ + 
2A^') - 2AAj/ ) + 1 or + 2A0& / ) - 2AA < g/ ) + 1. 

VI. CONCLUSION 

The problem of two-dimensional turbulence was investigated, taking into account the 
presence of three-dimensional perturbations. They were introduced in an effective way, 
represented by random forcing terms which act at small scales in the two-dimensional Navier- 
Stokes equations, as suggested by experimental observations. A coupling constant, g, related 
to the strength of these additional forces, allowed us to write a power expansion for the 
velocity field, containing also a compressible part. An infinite set of equations was found by 
just collecting terms with the same powers of g. The components ip^ an d <p^ p \ appearing 
in the power expansion of the velocity field were assumed to be primary operators of some 
conformal minimal model. We obtained, then, from point-splitted products of operators, 
a group of conditions in order to have a solution of the Hopf equations and to reproduce 
the situation of a constant enstrophy or energy flux through the inertial range. In the 
constant flux conditions, large and small scale terms were defined and evaluated by means 
of an extension of the dimensional argument employed formely in the study of analogous 
correlation functions. An analysis of the first six minimal models of unperturbed conformal 
turbulence was performed, showing that the picture of a constant enstrophy cascade is in 
good agreement with experimental data, yielding inertial range exponents, for the strong 
coupling regime, g » 1, very close to the ones observed in the laboratory. Regarding the 
energy cascade case, we noticed that most of the minimal models considered in our study 
were unable to give solutions for the perturbed system. Only one solution was obtained, 
with inertial range exponents around —3.0. It would be interesting to investigate further 
minimal models for the energy case, in order to see if a closer connection with the results 
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indicated in experiments could be reached. 

From tables 1 and 2, and figure 1, we see clearly that there are many solutions, differing 
by just one of the fields ip\ or <po, which give exactly the same inertial range exponents, in the 
strong coupling regime. It is tempting, then, to conjecture that one could find "plateaux" 
for the spectral slopes, while varying some set of external parameters. This question is 
contained, of course, in the deeper problem of how to match large scale properties of the 
fluid with the minimal models describing the inertial range. 

A point which deserves attention is the crossover between unperturbed conformal tur- 
bulence and the results obtained in the strong coupling regime. A bridge between these 
two situations may be investigated not only by varying g, as we did in section IV, but 
also through fi — > 0, when the effects of small scale three-dimensional perturbations on the 
constant enstrophy or energy fluxes become negligible. Finally, it is important to stress 
that a standard direct numerical simulation of equations ( |3.4| )- (|3ij| ), up to some level n in 
their hierarchy, would be an interesting way to study the above questions and the physical 
assumptions addressed in the present work. 
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TABLE AND FIGURE CAPTIONS 



Table 1. Solutions for the constant enstrophy flux condition. The first six models of 
unperturbed conformal turbulence were analyzed, all of them yielding possible definitions 
of ipi and 0o- Here we show the results for the minimal models (2,21), (3,25) and (3,26). 

Table 2. Solutions for the constant energy flux condition. The analysis of the first six 
models of unperturbed conformal turbulence showed that most of them were "blocked" by 
the presence of perturbations. The only obtained solution corresponds to the model (10,59). 

Table 3. Statistical data related to the solutions found for the constant enstrophy flux 
condition, in the strong coupling regime, where a comparison is made with the unperturbed 
values of the inertial range exponents. 

Figure 1. Graphic representation of the inertial range exponents in the enstrophy cas- 
cade case. Figures a), b) and c) refer to the minimal models (6,55), (7,62) and (8,67), 
respectively. The horizontal axis, labelled by (m,n), represents the most relevant field be- 
tween ipi and 0o, in the strong coupling regime. The ordering of fields is the same as in the 
tables. 
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Table 1. 



minimal model: (2,21) 




00 


4A^i + 1 


4A0o + 1 




00 


4A^i + 1 


4A0 O + 1 


(1,10) 


(1,6) 


-5.12 


-4.8 


(1,9) 


(1,9) 


-7.38 


-7.38 


(1,10) 


(1,7) 


-5.12 


-5.24 


(1,9) 


(1,10) 


-7.38 


-7.57 


(1,10) 


(1,8) 


-5.12 


-5.44 


minimal model: (3,25) 


(1,10) 


(1,9) 


-5.12 


-5.4 




00 


4A^! + 1 


4A0 O + 1 


(1,10) 


(1,10) 


-5.12 


-5.12 


(1,6) 


(1,6) 


-4.8 


-4.8 


(1,10) 


(1,H) 


-5.12 


-4.6 


(1,6) 


(1,7) 


-4.8 


-5.24 


(1,11) 


(1,6) 


-4.6 


-4.8 


(1,6) 


(1,8) 


-4.8 


-5.44 


(1,11) 


(1,7) 


-4.6 


-5.24 


(1,6) 


(1,9) 


-4.8 


-5.4 


(1,11) 


(1,8) 


-4.6 


-5.44 


(1,6) 


(1,10) 


-4.8 


-5.12 


(1,11) 


(1,9) 


-4.6 


-5.4 


(1,6) 


(1,H) 


-4.8 


-4.6 


(1,11) 


(1,10) 


-4.6 


-5.12 


(1,7) 


(1,6) 


-5.24 


-4.8 


(1,11) 


(1,H) 


-4.6 


-4.6 


(1,7) 


(1,7) 


-5.24 


-5.24 


minimal model: (3,26) 


(1,7) 


(1,8) 


-5.24 


-5.44 




00 


4A^i + 1 


4A0o + 1 


(1,7) 


(1,9) 


-5.24 


-5.4 


(1,6) 


(1,6) 


-4.96 


-4.96 


(1,7) 


(1,10) 


-5.24 


-5.12 


(1,6) 


(1,7) 


-4.96 


-5.46 


(1,7) 


(1,H) 


-5.24 


-4.6 


(1,6) 


(1,8) 


-4.96 


-5.73 


(1,8) 


(1,7) 


-5.44 


-5.24 


(1,6) 


(1,9) 


-4.96 


-5.77 


(1,8) 


(1,9) 


-5.44 


-5.4 


(1,6) 


(1,10) 


-4.96 


-5.58 


(1,8) 


(1,11) 


-5.44 


-4.6 


(1,6) 


(1,H) 


-4.96 


-5.15 


(1,9) 


(1,6) 


-5.4 


-4.8 


(1,7) 


(1,6) 


-5.46 


-4.96 


(1,9) 


(1,8) 


-5.4 


-5.44 


(1,7) 


(1,7) 


-5.46 


-5.46 


(1,9) 


(1,10) 


-5.4 


-5.12 


(1,7) 


(1,8) 


-5.46 


-5.73 
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00 


4A^i + 1 


4A0o + 1 


(1,7) 


(1,9) 


-5.46 


-5.77 


(1,7) 


(1,10) 


-5.46 


-5.58 


(1,7) 


(1,H) 


-5.46 


-5.15 


(1,8) 


(1,7) 


-5.73 


-5.46 


(1,8) 


(1,9) 


-5.73 


-5.77 


(1,8) 


(1,H) 


-5.73 


-5.15 


(1,9) 


(1,6) 


-5.77 


-4.96 


(1,9) 


(1,8) 


-5.77 


-5.73 


(1,9) 


(1,10) 


-5.77 


-5.58 


(1,10) 


(1,6) 


-5.58 


-4.96 


(1,10) 


(1,7) 


-5.58 


-5.46 


(1,10) 


(1,8) 


-5.58 


-5.73 


(1,10) 


(1,9) 


-5.58 


-5.77 


(1,10) 


(1,10) 


-5.58 


-5.58 


(1,10) 


(1,H) 


-5.58 


-5.15 


(1,11) 


(1,6) 


-5.15 


-4.96 


(1,11) 


(1,7) 


-5.15 


-5.46 


(1,11) 


(1,8) 


-5.15 


-5.73 


(1,11) 


(1,9) 


-5.15 


-5.77 


(1,11) 


(1,10) 


-5.15 


-5.58 


(1,11) 


(1,H) 


-5.15 


-5.15 
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Table 2. 



minimal model: (10,59) 




00 


4A^! + 1 


4A0o + 1 


(1,6) 


(1,6) 


-3.07 


-3.07 


(1,6) 


(2,12) 


-3.07 


-3.06 


(1,6) 


(3,18) 


-3.07 


-3.05 


(1,6) 


(4,24) 


-3.07 


-3.04 


(2,12) 


(1,6) 


-3.06 


-3.07 


(2,12) 


(2,12) 


-3.06 


-3.06 


(2,12) 


(3,18) 


-3.06 


-3.05 


(2,12) 


(4,24) 


-3.06 


-3.04 


(3,18) 


(1,6) 


-3.05 


-3.07 


(3,18) 


(2,12) 


-3.05 


-3.06 


(3,18) 


(3,18) 


-3.05 


-3.05 


(4,24) 


(1,6) 


-3.04 


-3.04 


(4,24) 


(2,12) 


-3.04 


-3.06 


(4,24) 


(4,24) 


-3.04 


-3.04 



Table 3. 



minimal model 


exponent (g=0) 


mean exponent (g ^ 0) 


standard deviation 


(3,25) 


-4.6 


-4.90 


0.28 


(3,26) 


-4.23 


-5.25 


0.27 


(6,55) 


-3.73 


-5.89 


0.21 


(7,62) 


-4.03 


-5.46 


0.28 


(8,67) 


-4.51 


-4.90 


0.34 
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